Special relativity in decays of hybrids 

Stanislaw D. Glazek 

AND 

Jakub Nar^bski 

Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, 

Poland 

(Received February 7, 2008) 

A decay of a heavy hybrid is expected to produce light mesons flying out 
with speeds comparable to the speed of light and phenomenological mod- 
els of the decay must respect symmetries of special relativity. We study 
consequences of this requirement in a class of simple constituent models 
with spin. Our models respect boost symmetry because they conform to 
the rules of a boost-invariant renormalization group procedure for effec- 
tive particles in light-front QCD. But rotational symmetry of the decay 
amplitude is not guaranteed and the parameters in the model wave func- 
tions must take special values in order to obtain the symmetry. When the 
effective interaction Hamiltonian responsible for a hybrid decay has the 
same structure as the gluon-quark-antiquark interaction term obtained by 
solving the renormalization group equations for Hamiltonians in first order 
perturbation theory, the non-relativistic image of a hybrid as built from a 
quark and an antiquark and a heavy gluon that typically resides between 
the quarks, cannot produce rotationally symmetric amplitude. However, 
there exists an alternative generic picture in the model that does satisfy 
the requirements of special relativity. Namely, the distance between the 
quark and antiquark must be much smaller than the distance between the 
gluon and the pair of quarks, as if a hybrid were similar to a gluonium in 
which one gluon is replaced by a quark-antiquark pair. 



PACS numbers: 13.25.Jx, 11.80.-m, 13.90. +i 



1. Introduction 

Special relativity symmetry imposes severe constraints on the constituent 
picture of decays of hybrids. When one attempts to construct a constituent 
model of hybrids based on the weak-coupling expansion in QCD £Q, treated 
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as a potential candidate to complement the strong-coupling lattice pic- 
ture 0E1; the parameters of the model cannot be treated as independent. In 
principle, QCD contains a small number of free parameters: quark masses, a 
coupling constant, and a scale parameter at which the masses and coupling 
are specified. In perturbation theory, the coupling depends on the scale 
through the ratio of that scale to Aqcd, the latter being adjusted in the 
scheme in which the constituent picture is being considered. However, not 
every scheme can clearly define the concept of constituents. It is particu- 
larly hard to describe constituent particles in the region of small virtualities 
where binding mechanism is at work. In fact, the constituent dynamics is so 
difficult to derive in QCD and solve precisely, especially for light quarks and 
gluons, that no clear quantum-mechanical picture in Minkowski space has 
been derived yet despite extensive studies. Thus, phenomenological images 
are based on models and such models contain additional parameters due to 
arbitrary simplifying assumptions that are not under control by a precise 
theory. 

Nevertheless, a phenomenological model constructed within a well-defined 
theoretical framework may help to finesse the leading approximation through 
agreement with data. Namely, if QCD is correct and does not require 
changes to precisely describe data, a model may hint at the structure that 
one should look for when attempting to solve the theory in a sequence of 
successive approximations. More importantly, data can provide constraints 
on models not only through discrete sets of numbers that correspond to the 
magnitudes of the considered quantities, but also through continuous sym- 
metries. The chief example is the Lorentz symmetry that includes boosts 
and rotations. The symmetry dictates the shape of functional dependence of 
a decay amplitude on the coordinates used to describe the outgoing particles. 
We show how this dictum works in a simple model. 

In the standard dynamical approach, rotational symmetry is kinematic 
(independent of interactions) and boost symmetry is dynamical (a change of 
frame of reference involves effects caused by interactions, such as a change 
in the number of constituents). Thus, in the standard approach, a model of 
a decay of a hybrid (or any other hadron) constructed in the center-of-mass 
frame of reference (CMF) of the hadron can easily respect the kinematical 
rotational symmetry in that frame of reference. However, it is not clear in the 
models that are based on the standard approach how well they respect the 
dynamical boost symmetry pQ. The boost symmetry would be very useful for 
checking the validity of all kinds of constituent models since QCD is supposed 
to yield only fully relativistic answers (HI El- Also, when one models a light 
meson in its own CMF as built from a fixed number of constituents and then 
uses the same constituent picture when the meson is moving with a speed 
close to the speed of light, which happens when a light meson is a product 
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in a decay of a hybrid, one has to verify if the model satisfies constraints 
that the boost symmetry imposes. But the constraints imposed by boost 
symmetry are hard to satisfy in constituent models in the standard dynamics 
because boost generators in the standard dynamics depend on interactions 
and change the number of constituents. 

In the light-front (LF) approach discussed here [Ej, boost symmetry is 
kinematical and it does not involve a change in the number of constituents 
no matter how fast a meson or a hybrid are moving. Instead, the rotational 
symmetry is dynamical. Therefore, if one wants to keep a fixed number of 
constituents, or when one seeks a physical picture in which the contributions 
from basis states with different than the leading number of constituents are 
small, it is the rotational symmetry of decay amplitudes (instead of boost 
symmetry in the standard approach) that begins to impose stringent con- 
straints on the LF wave functions that stand a chance to approximate so- 
lutions to QCD, if some constituent picture is actually valid in the theory 
Since the group of rotations is compact (the group of boosts is not) and it is 
already well-understood in non-relativistic quantum mechanics, the LF ap- 
proach provides an opportunity for studying dynamical constraints of special 
relativity in an intuitively familiar way. The rotational symmetry is a much 
more familiar concept of quantum theory than the boost symmetry. This 
is reflected in the fact that the latter symmetry is rarely discussed in terms 
of constituents in the context of QCD. But using the LF approach, one can 
gain some insight concerning structures that may emerge from a relativistic 
dynamical quantum theory by checking if a set of wave functions chosen in 
a model can provide a rotationally symmetric decay amplitude of a hybrid. 
We interpret our findings concerning hybrids in the LF approach from the 
point of view of models of non-hybrid mesons. 

Section discusses theoretical background for the present study. Our 
model assumptions are described in Section 03 Numerical results obtained 
from the symmetry constraints on the decay amplitude of a hybrid are de- 
scribed in Section Q] and discussed in Section 03 Section 03 contains our 
conclusions regarding the constituent structure of hybrids. 

2. Theoretical background 

A constituent model of a hybrid qualifies as a theoretically reasonable one 
when it is clear, at least in principle, how the constituents used in the model 
can be related to quarks and gluons in QCD within a single formulation of 
the theory. Of course, it is possible that QCD will be eventually solved and 
precise comparison with data will demonstrate that the theory requires some 
changes of currently unknown nature and implications [3J. Before this hap- 
pens, however, models of hybrid decays can be regarded as reasonable from 
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the theoretical point of view if they are designed in agreement with some 
framework for solving QCD. But since exact solutions are missing and the 
constituent picture for hadrons continues to pose major conceptual problems 
(there is no dynamical explanation in QCD of the phenomenological success 
of the constituent quark model, examples of studies can be found in 
there exists today a considerable room for varying parameters and adjust- 
ing them to data even in reasonable models. This status of theoretically 
reasonable models of hadrons will continue until QCD is solved precisely 
enough to derive the model pictures and remove the freedom in choosing 
their parameters. One should also keep in mind that a model may be suit- 
able for representing a theory for one type of data and not so reasonable 
for another type. There is an exception to this ambiguity: all theoretically 
reasonable models of hadrons must respect symmetries of special relativity. 
The alternative is to explain why the constraints of that symmetry can be 
ignored. 

The requirement of Lorentz symmetry imposes constraints on the wave 
function of the hybrid as soon as one has a candidate for the dynamical 
mechanism of the decay. This mechanism should also be derived from QCD. 
Unfortunately, constituent dynamics is not understood yet in terms of the 
theory. In the LF QCD, one can apply the renormalization group procedure 
for effective particles (RGPEP, see e.g. which shows the method in a 
considerably simpler case of heavy quarkonia, rather than the very difficult 
case of hybrids) and derive interactions of constituent quarks and gluons 
order by order in perturbation theory assuming an extremely small coupling 
constant, as if Aqcd were much smaller than it actually is. Certainly, when 
Aqcd is set to a realistic value, and when the appropriate renormalization 
group parameter is lowered down to the scale of the binding mechanism, 
the asymptotically free coupling constant increases to values for which the 
formal perturbative expansion in powers of the coupling constant is not 
expected to work well. But numerical tests in asymptotically free models 
with bound states suggest that perturbation theory in RGPEP may be 
able to identify parts of the structure of the interaction terms in effective 
Hamiltonians that dominate in the bound-state dynamics, using as small a 
coupling constant as one wishes in the process. When a structure is already 
found, the main effect is the increase of the coupling constant in front of the 
identified structure ^2]. This is true provided that one does not lower the 
renormalization group parameter in RGPEP too much and the calculation 
of the effective interaction is not cutting into the mechanism of formation of 
the bound states of interest. 

The first term relevant to hybrid decay that one derives in perturbation 
theory is the term in which an effective gluon turns into a pair of an effec- 
tive quark and an effective antiquark. Therefore, we assume here that the 
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interaction that leads to the decay of a hybrid is that of a constituent gluon 
decaying into a pair of a constituent quark and antiquark. The constituents 
are identified with effective particles in RGPEP. According to RGPEP, the 
relevant structure of the interaction is exactly like in the canonical Hamil- 
tonian for LF QCD except for two features. One feature is that the annihi- 
lation operator for the effective gluon, the creation operator for the effective 
quark, and the creation operator for the effective antiquark, all correspond to 
a small renormalization group scale A just above the scale where the binding 
mechanism is active. Another feature is that the interaction vertex contains 
a form factor f\, instead of being local as in the canonical theory that one 
starts with. The effective interaction term is denoted by Hj\. 

Once the interaction term Hj\ responsible for the decay of a hybrid 
state, denoted by \h), into two mesons, denoted by \p) and \b) (the letter p 
is chosen for a light meson, like meson n, and letter b is chosen for a much 
heavier meson, like meson b\ of mass 1235 MeV), is specified, the decay 
amplitude, denoted by A, is evaluated using the formula A = (pb\Hi\\h) . 
The symmetries of the decay amplitude A depend on the shapes of the wave 
functions of the hybrid, meson p, and meson b. Our question is: For what 
wave functions of the hybrid and two mesons one can obtain a spherically 
symmetric decay amplitude A for a ++ hybrid using the interaction term 

We choose the wave functions to correspond to the phenomenological 
images that underlie constituent models of hadrons. For example, we choose 
a Gaussian wave function of the relative momentum of quarks to model a 
meson. We demand that the width of that function is on the order of masses 
of the involved particles. Similar Gaussian wave functions are introduced for 
the hybrid state. There are also spin dependent factors for quarks and gluons 
in the states we consider that were not studied before pj. Our study includes 
several choices of these factors. The factors are built from the spinors and 
Dirac matrices which appear in the current operators that formally can 
produce meson or hybrid states with the quantum numbers we consider. We 
check the decay amplitudes of ++ hybrids into two types of mesons: two 
scalar or two pseudoscalar ones. The interaction term we use here differs 
from the scalar-gluon term used in the previous study jSj by inclusion of the 
gluon spin as dictated by QCD. Typically, the resulting decay amplitudes are 
not spherically symmetric. However, the degree of violation of the spherical 
symmetry depends on the values of the parameters we introduce in the wave 
functions and the RGPEP parameter A. But one can vary the parameters 
and check if there exist any choices for which an amplitude is spherically 
symmetric. It turns out that such choices do exist and we find them here 
by minimizing the deviation of A from spherical symmetry. 

The main assumption that is tested in our study is that the number of 
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constituents can be minimal. We know that the effective particle dynamics 
derived using RGPEP in LF QCD includes Hamiltonian terms that change 
the number of effective quarks and gluons. In the standard formulation 
of particle dynamics that evolves in time t = x° instead of the LF x + = 
x° + x 3 , special relativity requires that states of hadrons are built from Fock 
sectors with different numbers of virtual quanta. Even the vacuum state, 
with no hadrons at all, appears to be a very complex state whose structure 
eludes efforts of physicists to explain it. But the situation is different in 
the LF dynamics. When QCD is regulated in transverse (perpendicular to 
the z-axis) and longitudinal (along the front) direction, there is no creation 
of quanta from the bare vacuum and hadrons can be considered using an 
expansion into their Fock components. Moreover, the effective interactions 
that are obtained from RGPEP contain the vertex form factors f\ that may 
have a small width A in momentum space. These form factors prevent the 
interaction terms from easily producing additional constituents, even if the 
coupling constant is not small in comparison to 1. But since we do not 
know if one can approximate the solution for hadronic states in LF QCD by 
keeping only the smallest possible number of constituents with some A, the 
critical question is if there exists any reasonable choice of the parameters in 
a model with only a minimal number of constituents for which the model 
renders a spherically symmetric A. If such choices exist, what do we learn 
about the allowed model parameters from the symmetry requirement? 

It will be shown that all the sets of parameters that we obtain share some 
features. Some of these features turn out to be independent of all details in 
our treatment of spin of the effective quarks and the gluon. The conclusion 
obtained earlier in Ref. p], using a model with a spinless "gluon," is shown 
to be also valid when one includes spin. We study several options to do 
so and we find that the tendency observed in for a spinless gluon is a 
generic phenomenon, even though some changes do occur. But the general 
conclusion is that the probability distribution for the constituent quarks and 
a gluon in a scalar hybrid must resemble a state built from the gluon and 
an octet diquark. The diquark has a smaller size than the typical distance 
between the gluon and the diquark. The structure can be imagined as a 
gluonium with one gluon replaced by a small quark-antiquark pair. This is 
a result of pure fitting of the model parameters. Our study does not answer 
the question if or how this picture may arise from the effective LF dynamics 
in QCD. 

Although it is known how the effective dynamics can be derived order by 
order from QCD, the resulting eigenvalue equations for mesons or hybrids 
are too complex and the number of basis states too large for solving the 
equations completely without some guiding rules for simplifying the math- 
ematics. What we observe here is that the constraints of special relativity 
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strongly limit the acceptable wave functions if one assumes that physical 
states are dominated by the Fock sectors with the smallest possible numbers 
of effective constituents. The constraints of relativity force the wave func- 
tion parameters to take values that suggest a dominant role of gluons in the 
distribution of matter inside hybrids. Gluons seem to dictate to quarks what 
the latter must do, rather than vice versa, i.e., not as constituent models 
based on the picture for glueless hadrons suggest. Our main point is not 
that our model must be correct, but that the constraints of relativity on LF 
models with a minimal number of constituents are quite restrictive, can be 
implemented in practice, and point in new directions. 

Let us add that the problems with Lorentz symmetry in constituent 
models of bound states of quarks and gluons occur not only when an out- 
going meson is light and has to move fast in the rest frame of a decaying 
hybrid. They also occur when one considers a decay of a hybrid in fast mo- 
tion, which happens whenever the decay is a part of a bigger process that 
includes production and propagation with a high speed of the hybrid itself. 
Such circumstances may be of interest, for example, in a photoproduction 
of hybrids in motion. 



3. Assumptions 

The model we discuss is an extension of the scalar model from Ref. jH] 
and we adopt notation used there without changes. The new element here 
is the spin of a constituent gluon. In Ref. (SJ, gluons were treated as scalar 
particles. Here, the interaction Hamiltonian Hj\ that is responsible for the 
decay of the constituent gluon is taken directly from LF QCD with the 
RGPEP width parameter A near the scale of hadronic masses: 

n IX = gf x ij x ^A a ^ x . (1) 

We display below details of the term that creates a pair of an effective quark 
and an effective antiquark from an effective gluon. This is the only term 
that counts in our calculation of the decay amplitude A. The interaction 
term contains the vertex form factor fx- If we denote the invariant mass of 
the quark-antiquark pair by M. q q and the gluon mass by M g , then 

3.1. qq Mesons 

The qq meson wave functions are of the same type as in Ref. p]. 

\p) = J2 A 12 ! P + «(l + 2-p)^ PO (l,2)6l 1 4|0) , (3) 

12 J 
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where ^ PC (1, 2) is a product of color, flavor (isospin), spin, and momentum 
dependent factors: 

^ PC (1,2) = xlCpXc, xll P X l2 x\M l ^)Xs 2 ^p(l,2) , (4) 

with C p = l/v3 (color singlet), I p = l/v2 (isospin singlet), and S p (l,2) is 
a 2 x 2 spin matrix, sandwiched between two-component spinors. 5 denotes 
167r 3 times a three-dimensional <5-function of plus and transverse momenta 
of the particles indicated in the argument, S(k) = 16n 3 5(k + )d^ (fe -1 ), see [E|. 
The wave function ip p (l,2) is chosen to be Gaussian function, 

tp p (l,2) = N p N pm (k 12 ) N ps (k 12 ) exp 

where k\ 2 is the relative three-momentum of the quarks in their center of 
mass system, see AppendixEJ The additional functions N pm and N ps p] are 
introduced entirely ad hoc. One option we investigate is that these functions 
are kept equal 1. In this case, the momentum integrals involve the relativistic 
momentum-space measure and full complexity of factors resulting from the 
relativistic spin structure. In particular, the normalization of a meson state 
is given by an integral of a function that is a product of the square of 
the Gaussian function, a factor resulting from the relativistic momentum- 
space measure, and a complex momentum-dependent spin factor. The other 
options we investigate are that the functions N pm or N ps are chosen to 
cancel the relativistic momentum-space measure or the spin factor in the 
normalization integral, respectively. The normalization condition is {p\p') = 
p + 8(p — p'). When both the measure and spin factors are canceled by N pm 
(measure) and N ps (spin), the normalization integral is a plain Gaussian 
integral as in a non-relativistic quantum mechanics. We investigate these 
options to find out how strongly the relativity constraints on the amplitude 
A depend on different factors. The same type of factors as N pm or N ps in 
the meson p are introduced in the meson b and denoted by or iV& s . 
All factors N are listed in the Appendix EI The cases we discuss here are 
described as N = 1 (full relativistic complexity of the model wave functions) 
or N ^ 1 (non-relativistic appearance of the normalization integrals of the 
model wave functions). 

For J PC = ++ mesons, we have the following spin factor 

xt 1 Sp(l,2)x S2 =u 1 v 2 , (6) 

(the notation is the same as in Ref. j^]), where u\ and v 2 are Dirac spinors 
for quarks. For J PC = ^ mesons (pseudoscalar mesons) we use 

^^(1,2)^=^. (7) 



-k 2 
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The LF spinors we use here are described in the Appendix El 

The above model meson states require comments concerning their pa- 
rameters, spins, and quark content. First of all, we use the name 6-meson 
here to indicate generically that the meson is relatively heavy, like mesons 
bi or 77(1295) or others in the same range of masses. The spin of the meson 
b is assumed to be zero, and we consider only scalars and pseudoscalars to 
find out the consequences of the special relativity constraints in simplest 
and most transparent cases. Thus, we do not describe the spin of real b\ 
mesons. The parameters of the wave function of a 6-meson in our model are 
not constrained to explain properties of any real meson and they are left free 
within a considerable range in order to check what, if any, combination of 
all parameters can produce rotationally symmetric decay amplitudes. The 
issue here is not if we can fit a model to data when we ignore gluons in a 
model of ordinary mesons, but if ignoring gluons in a potentially valid ef- 
fective constituent picture in QCD is allowed by special relativity symmetry 
even in principle, before a dynamical analysis is attempted. 

The same applies in the case of a light meson, called here p-meson, which 
can have a mass as small as a 7r-meson. However, in the case of a pion, the 
assumption that such light meson is dominated by a quark-antiquark compo- 
nent may be considered merely a mock up when one attempts to understand 
the structure of light mesons in terms of canonical (almost massless) quark 
degrees of freedom in QCD, where the problem of breakdown of chiral sym- 
metry requires a careful statement, or when one tries to create a strong 
binding effect in a naive potential model. We are facing a problem that on 
the one hand hadrons can be classified in terms of constituent quarks and 
lightest mesons belong to the same scheme as the heavier ones, and on the 
other the phenomenon of chiral symmetry breaking in canonical QCD is not 
explained quantitatively. Our model study of symmetry constraints in the 
effective constituent picture based on RGPEP in QCD is not solving this 
problem. What matters is that the effective particle picture is not necessar- 
ily wrong. Namely, the effective quarks may have large masses and chiral 
symmetry may be already explicitly broken in the effective Hamiltonian that 
has the width A comparable with hadronic masses. At the same time, the 
binding potentials in H\ in QCD may differ in the case of 7r-mesons and 
in the case of heavier ones. The key question here is not what dynamical 
mechanism might be responsible for the success of the constituent classifica- 
tion of hadrons, whether it is vacuum condensates in standard dynamics or 
corresponding special terms in a LF Hamiltonian, but if such minimal con- 
stituent picture can satisfy constraints of special relativity, assuming a most 
plausible Hamiltonian term that can lead to the decay of a hybrid. There- 
fore, we take the stance at this stage of the development that anything goes 
that produces relativity with constituents and we ask if this condition can 
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be satisfied in any, even remotely plausible way from the point of view of 
currently popular models. Our major issue is if the LF QCD constituent 
picture has a chance to obey rotational symmetry. We do not solve the 
dynamics of formation of pions here and we will call the light meson p, not 

7T. 

3.2. Hybrid Meson 
Our model for ++ hybrid state is of the form (cf. p]) 

\h) = I [123] h+~S(l + 2 + 3 - h)^ PC (l, 2, 3) 6^44310) . (8) 

123 J 

The creation operator for an effective gluon of width A, a^ 3 , carries spin, and 
the wave function ^ h jPC {\, 2, 3) depends on this spin. The wave function is a 
product of color, flavor (isospin), spin, and Gaussian functions of the relative 
momenta of the three particles: 

* JTO (1,2,3) = xlC c h a Xc 2 xljhX l2 xl 1 ^(l,2,3) Xs2 V/ l (l,2,3) , (9) 

with C^ 3 =t C3 /2 and Ih = l/\/2- The momentum dependent factor (hybrid 
wave function) is assumed to have the gg-cluster form EH 

^(1,2,3) =N h N hm (kg,k g ) N hs (k q ,k g ) exp 

with typical Gaussian functions of the relative momenta. Namely, k q is 
the three-momentum of the quark in the CMF of the quark-antiquark pair, 
and kg is the three-momentum of the gluon in the CMF of the quark, an- 
tiquark, and gluon (see Appendix El for details). Since we use the LF form 
of dynamics, the CMFs are well defined and separation of the relative and 
CMF motion for any state is purely kinematical. The optional additional 
factors Nhm and Nh s can be again kept equal 1 or chosen so that the wave 
function normalization condition has a non-relativistic appearance for three 
constituents similarly to the wave functions of mesons, see AppendixEl For 
consistency, we will always either introduce all factors ./V for both mesons 
and hybrid equal 1 (the cases labeled N = 1 ), or insert all factors iV such 
that all our states are normalized through the same integrals as in a non- 
relativistic theory (the cases labeled N ^= 1). 

The simplest scalar-hybrid spin factor that may be considered in a rela- 
tivistic theory is 

XiSh(h 2, 3)x2 = ui-f^v 2 , (11) 



kg 




kg 


w 2 hq 


exp 


wl 9 
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where £3 is the gluon polarization four vector. This case will be referred to 
as the case u$v. For the effective gluon in the light-front gauge A + = 0, the 
polarization vector has the following components: 



The sum over gluon polarizations is 



E e tA = sT + 3 ,+ 3 • (13) 

In the above formula, the component k% of the gluon momentum is the same 
as for massless gluons, = k^ 2 /k^ . But in evaluating kinetic energy of 
the effective gluons, we introduce the gluon effective mass parameter m g 
that can depend on the RGPEP parameter A. We do not know the value of 
m g and we leave it as a free parameter in our Gaussian wave function. 

An alternative structure for the spin factor, pertaining to non-abelian 
gauge symmetry but not necessarily better than Eq. Ill 111 from the dynamical 
point of view in the effective theory, is 

xl S h (l, 2, 3) Xs2 = ui^v 2 G» V P V , (14) 

where G^ u = k^e v ka — k^e ka and P = k\ + A; 2 + A; 3 . The four- vector e v ka is the 
polarization vector for massless gauge bosons. But k% can be calculated as 
if the gluon mass were 0, or using the parameter m g , and we do not know 
which way is more realistic in a dynamical theory. Therefore, we insert the 
unknown mass m g in the formula k^ = (k% 2 + m 2 g )/k^ and check for what 
values of m g the resulting decay amplitude A is spherically symmetric in 
the hybrid CMF. This case is referred to as uGPv. 

We also consider alternative versions of the spin factor, where we calcu- 
late k^ as if the gluon mass should be kept in G^ v and/or in P". These 
cases are referred to as uGPv, uGPv, and uGPv, respectively. The tilde 
means that we put m g = in evaluating k^ in the factor that is labeled 
with the tilde. 

The alternative spin factors resemble the one that occurs in the operator 
structures used in lattice calculations ^HE]- In the hybrid CMF, where 
P = 0, this factor reduces to u^vGio, which is is a combination of the 
components of the quark current and the chromoelectric gluon field. 

Note that the gluon momentum component k% = (k^ 2 + m 2 )/k^ does 
not contribute to the Lorentz product k^e^ in the chosen gauge and the en- 
tire small group of the Poincare transformations that preserve the light-front 
hyperplane also does not change the gauge condition A + = 0. Thus, we can 
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safely boost the hybrid state using kinematical relations and our assignment 
of mass to the effective gluon does not interfere with our choice of gauge. 
This is important because our model would not be reasonable otherwise. 
Namely, if the small group and gauge choice would not commute, we would 
not be able to construct the states of mesons and hybrids in motion without 
changing the gauge. The latter change would be associated with altering 
interaction terms in the effective Hamiltonian of width A. Fortunately, our 
choice of the LF dynamics, instead of the standard one, offers a possibility 
of keeping boost invariance using one and the same choice of gauge in all 
frames of reference that can be reached by the boosts. Note also that our 
RGPEP procedure respects this commutativity because it is invariant under 
the small group and our model is reasonable in this respect. 

As explained earlier, in contrast to the standard approaches where ro- 
tational symmetry is kinematical and boosts are dynamical, in the LF ap- 
proach the rotational symmetry is dynamical. Therefore, we now have to 
check to what extent our models can guarantee that the resulting decay am- 
plitude is spherically symmetric in the rest frame of the hybrid. We do this 
in the next section using both Eqs. ffTTf) and lfl"4")l in a number of cases that 
we have introduced above. 

4. Symmetry constraints 

The decay amplitude of a scalar hybrid ( J PC = ++ ) into two mesons, 
either two J PC = ++ mesons or two h mesons, should be spherically 
symmetric. But a constituent model built in the LF scheme introduces 
dependence on the angle 6 between the z-axis and the direction of flight of 
the light meson. This effect is the price we pay for boost invariance. The 
effect was discovered and initially studied in a scalar model in Ref. p|. 

Fig. H illustrates how badly the rotational symmetry is violated in a 
decay into two scalar mesons (spin factors uv) when the light meson mass 
varies from 664 MeV toward the value of 138 MeV. In this figure the factors 
JV), m , Nh s , Nb m , Nb s , N pm , and N ps , in the wave functions, are kept different 
from 1 to secure non-relativistic normalization. In the case of 664 MeV, the 
hybrid mass rah is just above the threshold of m& + m p and the product 
mesons can barely move. The LF constituent model renders an amplitude 
that does not depend on the angle 6. In the 138 MeV case, corresponding 
to the 7r-mesons, the hybrid mass is far above the threshold and the light 
meson has a highly relativistic velocity. In this case, the amplitude depends 
on the angle 6 to an unacceptable degree. The effect is caused by the fact 
that when the outgoing light meson flies against the z-axis with nearly speed 
of light, it must be built from quarks that have very small momentum k + 
and such quarks are suppressed in the case of wave functions used in the 
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Fig. 1. The angular dependence of the decay amplitude for various masses of the 
light meson, m p . The plot a) is for the hybrid spin factor equal uj^ve^, and b) 
is for spin factor uj^vG^P^- In both cases the hybrid meson decays into two 
J PC = ++ scalar mesons (spin factors uv). The model wave functions contain 
factors N ps , N pm , Nb s , Nb m and Nh s , Nh m that secure that the normalization 
integrals have a non-relativistic appearance of integrals of plain Gaussian functions 
(case TV 7^ 1). All parameters of the wave functions are given in the first column 
of Table □ 
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Fig. 2. The same as in Fig.QJ but with the factors N = 1. 



calculation. The parameters that we use are given in the first column of 
Table H on page ED (cf. 0)- Figs. [21 to Q] show the same effect, but with 
all factors Nh mi iV^ s , A^ m , iV^, N pm , and N ps equal 1, or in the case of 
two pseudoscalar (J PC = h ) mesons instead of the scalar ones, or in the 
case where the two changes are combined. In order to satisfy constraints of 
special relativity, the decay amplitude should not depend on the angle 0. 
Broadly speaking, the violation of rotational symmetry results from the 
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Fig. 3. The same as in Fig. but for hybrid decay into two J PC = + pseudo- 
scalar mesons (spin factors u"f 5 v). 
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Fig. 4. The same as in Fig. but with the combined effect due to the changes 
from Fig. |2 and El N = 1 and the hybrid decays into two pseudo-scalar mesons. 



fact that the model wave functions are not constrained dynamically by any 
underlying relativistic theory. Given that the model is reasonable, in the 
sense that (1) the meson states are formed using well-defined degrees of free- 
dom that appear in the LF Hamiltonian H\ in QCD with a small RGPEP 
parameter A, (2) the boost symmetry is preserved exactly, and (3) the decay 
is driven by an interaction term in the same Hamiltonian, the most ques- 
tionable element of the model is the assumption that a small number of con- 
stituents is sufficient to build a solution of a relativistic theory. The model 
assumes that the number of constituents is the smallest possible. It may fail 
to produce rotational symmetry because the symmetry is dynamical in the 
LF scheme and the interactions can change the number of constituents 0[Bj- 
We see in Figs. ^ to 0] that all models we test respect rotational symmetry 
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very well in non-relativistic decays. But in the relativistic decays, all the 
models fail more or less equally badly (on average, decays into pseudoscalars 
are a bit less wrong than decays into scalars because pseudoscalars are domi- 
nated by non-relativistic, momentum-independent components in their wave 
functions). Does this mean that all models based on the assumption of the 
smallest number of constituents must be entirely wrong? 

We find that the answer to this question is no: a minimal constituent 
model does not have to be wrong. Since the rotational symmetry is dy- 
namical in the approach we study, it is not known if the parameters listed 
in Table ^ in the first column do correspond to a solution of a relativistic 
theory. Suppose that a different set of parameters should be used in a rea- 
sonable model that approximates a solution of a relativistic theory. Can 
one find a set of parameters in the constituent wave functions for which 
the required rotational symmetry of the decay amplitude is obtained? This 
question is found to have a positive answer but the sets of parameters that 
we find point to a new picture for the hybrids. The picture seems to be 
generic in the sense that its dominant features are independent of how the 
spin of the gluon and the spins of quarks are treated. Our numerical studies 
produce examples of models with a smallest number of constituents in which 
the rotational symmetry is respected well when one allows the parameters in 
the wave functions and the RGPEP scale A in the Hamiltonian to vary. The 
reader should remember that the number of variable parameters in the class 
of models we consider is 7 and there exists a great number of possibilities 
to check, each demanding a multidimensional integration for every value of 
the angle 9. 

Fig. [HI shows how well the rotational symmetry can be restored by select- 
ing a different set of parameters in the wave functions (the sets corresponding 
to Fig. El are given in columns |EK-d in Table ^) . Curves "a" and "b" repre- 
sent the decay amplitudes for a hybrid meson with the spin factor given in 
Eq. (|TT1l (case u$v) and the wave-function parameters given in the second 
column of Table ^ Curves "c" and "d" on the same figure show the decay 
amplitudes for a hybrid with the spin factor given by Eq. (|THl (case uGPv) 
and with the parameters given in the third column of Table ^ Curves "a" 
and "c" are obtained with N ^ 1, and curves "b" and "d" are obtained with 
N = 1. 

The optimal choice of the parameters that one obtains from the condition 
of rotational symmetry includes f3 p about twice smaller than the typical 
value of 0.4 GeV in the first column of Table which corresponds to the 
size of a real meson ir (we will return to this issue below). We could also 
find other sets of parameters with even smaller (3 P and considerably smaller 
quark masses, a feature observed already in Ref. pj. The cases shown here 
are obtained by starting a Powell local minimization procedure ^1 EE3 
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Fig. 5. The hybrid decay amplitude as function of 9 in four cases: a) u^v, N ^ 1 
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uGPv, JV^ 1 and fourth column of Table d) uGPv, N = 1 and fifth column of 
Table n] Decays into two J PC = ++ mesons. 
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Table 1. The parameters of wave functions (in GeV) that we use in Figs. 11151 fin 
the same notation as in Ref. 0). s means that the parameter is the same as in the 
left neighboring column. Columns [IJi-d display results of local minimization using 
Powell's procedure ^1, starting from the values given in the last column (labeled 
"hs"), which is the same as the corresponding column in Table I in Ref. 0. 



starting from the hs (for "heavy-scalar.") set of parameters that was found 
in a scalar model in Ref. 0. The nomenclature refers to the relatively heavy 
scalar particles that played the role of quarks in Ref. p]- 

The cases we display here are characterized by quite good spherical sym- 
metry in comparison to other locally optimal choices which we were also 
able to identify but which displayed more variation with 9. In choosing the 
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minimization we also adopted a criterion that the resulting spherically sym- 
metric amplitude should not be many orders of magnitude smaller than in 
the case of the parameters in the first column of Table ^ Note, however, 
that the amplitudes in Fig.EJare, in fact, a whole order of magnitude smaller 
than in Figs.^^1 This indicates how important the constraints of relativity 
can be for analysis of data. We should also add that the size of the coupling 
constant in the interaction Hamiltonian that drives the decay was fixed as 
in Ref. and never changed in the fit. 



par am 


m q 


m g 


Pp 


A, 


Phg 


Phq 


A 


min 


0.1 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


max 


0.6 


2.0 


0.8 


1.6 


2.0 


8.0 


8.0 



Table 2. The limits on parameters of the wave functions (in GeV) that we imposed 
using the Adaptive Simulated Annealing algorithm [18| . 

It is clear that one should not consider an unbiased minimization of sym- 
metry violation as a most reasonable approach. The minimization should 
include additional constraints, including restrictions such as the radius of 
a meson p, and as much of the dynamical constraints as possible. But in 
order to impose correlations such as the ones coming from the radius, one 
has to be very careful about how one calculates the radius and if that cal- 
culation does obey requirements of special relativity, which is a problem in 
itself. Concerning the dynamical constraints, we were not able and not even 
interested in imposing any such constraints at this stage, because we were 
only searching for the answer to the question if any choice of the parameters 
could produce spherical symmetry, and it is interesting that even a minimal 
model can produce the symmetry of the quality as good as shown in Fig. 
We remind the reader that the gluon spin introduces functions of momenta 
that vary rapidly with angles and it was not clear at all that any choice of 
parameters could lead to a constant amplitude. But once it is established 
that such result is possible, one can make further observations based on a 
systematic search through the space of the parameters. 

The simplest and least restrictive way of setting bounds on the parame- 
ters of the models we test is to limit all of the parameters to fixed intervals 
around values that are considered reasonable. Such least restrictive pa- 
rameter bounds adopted in the minimizations described below are given in 
Table □ 

We performed a global minimization of departures from rotational sym- 
metry using different measures of how much a decay amplitude differs from a 
constant as a function of the angle 9: a standard deviation from the average 
value (sum of squares of deviations from the average value, labeled "stddev" 
in AppendixEI), or maximum of the modulus of the deviations from the av- 



18 



10" 1 




10 -5 I , 1 , 1 , 1 

60 120 180 

e [deg] 

Fig. 6. The hybrid decay amplitude as function of 9 for the decay into two J PC — 
++ mesons in different cases: a) u^v, N ^ 1, first column of Table b) u$v, 
N = 1, second column of Table c) uGPv, N ^ 1, third column of Table d) 
uGPv, N = 1, fourth column of Table H 



Fig. # 




Eb 


Efc 


Hi 


spin 






uGPv 


uGPv 


term 


N ^ 1 


N = 1 


N ^ 1 


N = 1 


m h 


1.9 


s 


s 


s 


m b 


1.235 


s 


s 


s 


m p 


0.1375 


s 


s 


s 


m q 


0.152 


0.21 


0.17 


0.155 


m g 


1.28 


1.07 


1.70 


1.90 




0.132 


0.219 


0.1 


0.1 


A, 


0.320 


0.536 


0.321 


0.371 


Phg 


0.766 


1.05 


0.267 


0.263 


Phq 


8.0 


7.83 


4.59 


5.73 


A 


7.68 


6.13 


2.86 


7.98 



Table 3. The optimal parameters of the wave functions (in GeV) for a decay 
into two J PC = ++ mesons. These are results of a global minimization using 
ASA pH], minimizing standard deviation from the average value of the amplitude, 
for parameters within limits in Table El The resulting amplitudes are shown in 
Fig. E| The bold face numbers are on the limit of the allowed range. 

erage value (labeled "maxdev" in the Appendix EJ, both measured relative 
to the average. Our global minimization within the assumed bounds is done 
using Adaptive Simulated Annealing (ASA) [18J. 

Fig. El shows results obtained using Eq. tfTTfl for the spin factor of a 
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Fig. 7. The hybrid decay amplitude as function of 9 for the decay into two J = 
++ mesons in different cases: a) uGPv with N ^ 1, b) uGPv with AT = 1, c) 
uGPw with 2V j= 1, d) uGPv with N = 1, e) wGPw with N ^ 1, i) uGPv with 
N = 1. Parameters are given in Tabled 
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Table 4. The optimal parameters of the wave functions (in GeV) for a decay 
into two J PC = 0++ mesons. These are results of global minimization using 
ASA |18j . minimizing standard deviation from the average value of the amplitude, 
for parameters within limits in Table [2j The resulting amplitudes are shown in 
Fig. [3 The bold face numbers are on the limit of the allowed range. 



hybrid meson: curve "a" in case N ^ 1 and "b" in case N = 1, and using 
Eq. ifTl^l : curve "c" in case N ^ 1 and "d" in case N = 1, all cases for a 
decay into two scalar (uv) mesons. The factors N have considerable impact 
on the magnitude of the amplitudes and can compensate or dramatically 
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Fig. 8. The hybrid decay amplitude as a function of 6, the same as in Fig. but 
for a decay into two J PC = h mesons (spin factor uj 5 v). The optimal wave 
function parameters are given in the correspondingly marked columns of Table El 



Fig. # 


Ek 


Eb 


Efc 


m 


spin 


ufv 




uGPv 


uGPv 


term 


N^l 


N = 1 


N + 1 


N = 1 


m h 


1.9 


s 


s 


s 


m b 


1.235 


s 


s 


s 


m p 


0.1375 


s 


s 


s 


m q 


0.15 


0.16 


0.18 


0.19 


m g 


1.08 


0.88 


1.69 


1.79 


Pp 


0.21 


0.25 


0.16 


0.22 


0b 


0.59 


0.59 


0.42 


0.30 


Phg 


0.59 


0.72 


0.68 


0.47 


Phq 


2.40 


2.62 


6.96 


6.80 


A 


3.97 


3.71 


2.29 


2.39 



The optimal parameters of the wave functions (in GeV) for a decay 

+ 



Table 5. 

into two J PC - o-+ mesons. These are results of global minimization using 
ASA [TS] . minimizing standard deviation from the average value of the amplitude, 
for parameters within limits in Table [2j The resulting amplitudes are shown in 
Fig. El 



enhance the effect of changing the spin factors. This means that one should 
probably not trust models of hybrids that are based solely on non-relativistic 
intuitions. 

Results for decays into two J PC = mesons (u^v spin factor) are 
shown in Fig. El The corresponding values of the wave function parameters 
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Fig. 9. The hybrid decay amplitude as function of 9 for the decay into two pseu- 
doscalar J PC = h mesons in different cases: a) uGPv with N ^ 1, b) uGPv 
with N = 1, c) uGPv with AT ^ 1, d) uGPv with AT = 1, e) uGPv with N 
uGPv with N = 1. Parameters are given in Table El 
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Table 6. The optimal parameters of the wave functions (in GeV) for a decay into 
two " + mesons. These are results of global minimization using ASA |18j . 

minimizing standard deviation from the mean amplitude for parameters within 
limits in Table The resulting amplitudes are shown in Fig. El 



are given in Table [3 Similarly, Figs. and show results for alternate 
choices of hybrid spin factor: cases referred as uGPv, uGPv and uGPv. 

In Ref. p], there were found two locally best sets of values of parameters 
in each of the two cases: one case with all constituents being scalars, and 
another one with fermionic quarks and a scalar gluon. These good sets were 
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Fig. 10. The hybrid decay amplitude, like in Fig. for the decay into two scalar 
mesons, as function of 6. Parameters are given in the corresponding columns of 
Tabled the lower limit on the quark mass is m q > 300 MeV. 



characterized by either light or heavy mass of the quarks (Table I and Fig. 5 
in Ref. pj). In our studies, including the spin of the gluon, the complete 
ASA algorithm finds only one best set of the parameters that minimizes 
deviation from rotational symmetry in every case we consider. These best 
sets appear with small quark masses and small j3 p . Such small j3 p implies a 
too large size of the meson p, apparently corresponding to a too weak binding 
of too light quarks, as if the number of such light effective constituents could 
not be only minimal. But we can find different minima when we impose an 
additional restriction that the quark mass, m q , is "heavy", i.e. greater than 
300 MeV. Other parameters are still limited to the intervals given in Tabled 

An example of such restricted minimization for "heavy" effective con- 
stituent quarks is shown in Fig. The corresponding optimal parameters 
are given in Table In all cases except the case "d," the size of the light 
meson p is now much closer to the size of the real mesons ir. This shows 
that the size (radius) of the light meson may not be as big an issue as one 
might think on the basis of a search for the best parameters allowing quarks 
to be much lighter than 300 MeV (this case was discussed earlier). 

Note that the assumption of case "d," that the spin factor in the hybrid 
wave function contains a four-momentum of a massive gluon, leads to a small 
optimal quark mass and a very small decay amplitude. The constraints of 
rotational symmetry promise to be very useful in future studies of reasonable 
models because when they are combined with inspection of observables such 
as radii one immediately obtains large differences between predictions based 
on different models. 
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Table 7. The optimal parameters of the wave functions (in GeV) from a global 
minimizaton using ASA ^Hl> minimizing standard deviation relative to the mean 
amplitude with all parameters within limits like in Table except for the lower 
limit on the quark mass, m q > 300 MeV. The resulting amplitudes are shown in 

Fig. ma 

5. Discussion 

First of all, let us note that the inclusion of spin of a constituent gluon 
does not change the previously obtained result for scalar "gluons" that 
rotational symmetry is restored when the quark-antiquark pair momentum- 
space width in the hybrid, (3h q , is about the same in size as the width A 
in the vertex form factor in the renormalized interaction Hamiltonian H\ 
in LF QCD, and both are on the order of 4-6 GeV, much larger than all 
other parameters. On the basis of our study of many cases with different 
ways of including the gluon spin and different cases of meson spin factors, 
we can state that the required relative momentum-space structure of the 
wave functions appears to be qualitatively independent of the spin of the 
effective constituents. Quite generally, the parameters f3h q and A appear to 
have to be about 5 times larger than the other parameters, see Tables ^ 
and Oh E This result suggests that the quark pair should be thought about 
as spatially small in comparison to the size of the hybrid. It also suggests 
that the pair could originate from a gluon that belonged to a gluonium before 
the interaction changed one gluon into the pair. 

One may worry that this result is obtained by minimizing just one ob- 
servable in a space of seven parameters. Here comes the strength of the 
continuous symmetry condition on a reasonable model: sepcial relativity 
provides infinitely many conditions instead of just one - the just one decay 
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amplitude must not be a function of the angle. It was highly questionable 
that it was possible to even come close to a constant function of the angle 
before we carried our study including the singular spin factors for gluons. 
The remarkable fact is that the symmetry cannot be satisfied with a mini- 
mal constituent picture unless the hybrid built from a pair and a constituent 
gluon looks differently than expected assuming that the gluons follow quarks 
and this picture does not seem to depend on any particular detail but only 
on the major assumption that a minimal constituent model can approximate 
the effective LF dynamics. 

Another worry concerns the small spatial size of the qq-p&ir, about 4 to 
5 times smaller than the distance between the gluon and the octet diquark, 
correlated with A on the order of 3 or even 5 GeV. The real question is for 
what values of A a constituent picture of hybrids may work. In principle, if 
RGPEP equations were solved exactly, no physical result should depend on 
A and no matter what A is used one should obtain rotationally symmetric 
decay amplitude provided that the decay is calculated exactly using exact 
solutions for the participating hadrons. But we know |l2j that in order to 
approximate the full dynamics of an asymptotically free theory by a simple 
picture one has to lower A to values that are about twice above the scale of 
eigenvalues one is seeking to describe. One cannot lower A to smaller values 
using perturbation theory in RGPEP because the resulting H\ would begin 
to contain too large errors due to cutting into the mechanism of binding. 
Thus, the scale we obtain from the heuristic fit is quite reasonable. On the 
other hand, one may worry that no tight diquark clusters are seen in the 
proton deep inelastic structure. One possible explanation of such special 
feature of hybrids could be that they contain octet diquarks interacting 
with constituent gluons, while in the proton we have primarily triplets and 
diquark antitriplet and no counterpart of the constituent gluon structure. 
Our symmetry study in hybrids should be seen as groping into the sectors 
of effective color dynamics that are squeezed out of and cannot be seen in 
nucleons. 

There exist small differences between the case of decays into two scalar 
and two pseudoscalar mesons, and between various choices for handling the 
spin of an effective gluon in the hybrid. But the outstanding feature that 
Phq is about the same as A and both are larger than the rest of parameters 
is common to all the cases we studied. Also, the mass of the gluon, m g , 
appears to have to be much greater than the masses of the quarks, m q . 

The hybrid spin wave function from Eq. (|T4"|) . inspired by the lattice 
operators but with a massive four-momentum for the gluon, gives minima 
with a much smaller decay amplitude A than the simplest hybrid spin wave 
function of Eq. (jTT|) . This result shows that one has to be very careful about 
treatment of spin of gluons in model building. 
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Let us stress that the absolute size of the amplitude is not under control 
in our calculation because we did not include the dynamical constraints 
between the coupling constant g, the RGPEP parameter A, and the wave 
function parameters. But there exists a systematic trend in all our results, 
which requires further study and is not understood here. Namely, when 
the parameters of the wave functions are varied from the values in the first 
column of Table to the values required by rotational symmetry, like in 
columns 2-5 of Table or the values in Tables OHEl and the size of the 
amplitude changes from about 0.9 to 0.1 in Figs. to the much smaller 
values of 10~ 2 or 10~ 3 in Fig.0 or 10" 2 or 10" 4 in Figs. ©to OH 

These are considerable changes in the order of magnitude. It seems 
unlikely that the coupling constant can vary by that much and compensate 
this change. The change is so large only because of the relativistic motion 
of the light outgoing meson. If the outgoing mesons are slow, rotational 
symmetry in our model is respected very accurately for typical values of 
parameters in non-relativistic constituent models, exemplified in the first 
column of Table ^ We are forced to conclude that a relativistic hybrid 
decay (including fast mesons p) may involve relativistic effects that are not 
accounted for in the non-relativistic phenomenology, cf. |19| . 

Another feature worth mentioning is that the minima for the wave func- 
tions motivated by lattice operator structures, Eq. (dJ, are narrower than 
in the case of Eq. ifTTjl . If the range of parameters for which a violation of 
rotational invariance is small is very narrow, the symmetry itself becomes a 
source of detailed information about the necessary values of the parameters 
even if the corresponding dynamical equations are too difficult to solve with 
comparable precision. By the same token, one obtains a very strict criterion 
for judgement of dynamical models that attempt to produce the relevant 
wave functions. 

6. Conclusion 

The example of a simple model described here shows that the wave func- 
tion parameters for hadrons involved in a relativistic decay of a hybrid must 
be strongly correlated in order that the decay amplitude satisfies require- 
ments of special relativity. Thanks to the use of the LF scheme, boost 
symmetry is respected exactly and the parameters are constrained by the 
condition of rotational symmetry. In the example, they have to take values 
that do not correspond to the picture based on the non-relativistic intuition 
that the gluons are mainly between two quarks. Instead, the relativistic 
effective constituent picture almost universally points toward the structure 
in which a heavy gluon is accompanied by a quark-antiquark pair that re- 
sembles a relatively small octet diquark. This is a stunning result because 
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it suggests that the picture with gluons playing a role of a relatively light 
chain, or a vibrating flux, or string between relatively heavy quarks may be 
not as realistic as one hopes for on the basis of non-relativistic intuition. 

The alternative hybrid structure occurs in a variety of cases that differ 
in details of the spin factors for gluons and quarks. But the requirement 
of relativistic symmetry turns out to be very restrictive when one demands 
that only sectors with the smallest possible number of constituents are im- 
portant. Therefore, we conclude that the effective constituent dynamics 
in QCD should be always considered including constraints of special rela- 
tivity. These constraints appear capable of forcing us to consider hadrons 
with significant gluon content not as if the gluons were just added to quarks 
and antiquarks, but as if gluons could actually dominante the dynamics of 
hybrids and force the quarks to adjust. 

Since the hybrid structure we are forced to seriously consider by the re- 
sults of this analysis contains a spatially tight octet diquark pair, as if the 
pair emerged from a constituent gluon through a single interaction in an 
effective LF QCD, one may ask if it is possible that such effective quark- 
antiquark-gluon states with a tight pair can mediate decays of usual mesons. 
A decay of a usual meson may proceed by an emission of a gluon from one 
quark and subsequent decay of the gluon into a new pair of quarks. The 
emerging two quarks and two antiquarks can form the mesons that are pro- 
duced in the decay of the usual meson. But if the pair accompanied by the 
intermediate gluon has to be small in size, as if the three effective particles 
had to form a structure similar to our finding for a hybrid, the intermedi- 
ate quark configuration would have to have a small overlap with the initial 
usual meson configuration. The decay mechanism through an intermediate 
hybrid meson would have small contribution to the total strong decay width. 
Would not this width be too small if the hybrid structure were as we obtain? 
Not necessarily, since in the effective theory there must exist other interac- 
tions that are capable of producing four effective quarks from two effective 
quarks. These interactions do not correspond to the intermediate excitation 
of a massive effective gluon and they are not characterized by the coupling 
of such gluons to quark-antiquark pairs. Examples of such interactions are 
present already in the canonical Hamiltonian in LF QCD. The canonical 
interactions are not mediated by emission or absorption of gluons, and they 
must contribute to the mechanism of strong decay of usual mesons in the 
effective theory characterized by width A on the order of hadronic masses. In 
addition, the RGPEP procedure generates more interactions that can turn a 
quark-antiquark pair into two such pairs without explicit creation and decay 
of a massive effective gluon corresponding to small A. Unfortunately, our 
study is not telling us anything about the dynamical structure of ordinary 
mesons and interactions that mediate their decays. It is limited to a pre- 
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liminary study of symmetry constraints in simplest models with a dominant 
hybrid component. 

The wave function parameters that we find to be preferred by the condi- 
tion of rotational symmetry of the decay amplitude of a model hybrid, may 
turn out to be invalid when the actual dynamics is included in the analysis. 
For example, it may turn out that the approximation by the Fock sectors 
with only the smallest possible number of constituents does not apply. But 
it is clear that the relativistic constraints cannot be ignored in the search 
for a leading constituent picture. 

Our discussion was limited to ++ hybrids for simplicity, while the most 
interesting from practical point of view is the structure of exotics ^3J 
EHII^ES]- One can change factors in the wave functions and change the 
quantum numbers of the hybrid states to the exotic values. For example, one 
can replace the color electric field by a color magnetic field, or introduce p- 
wave wave functions. Simple introduction of p-wave for gluon (introducing 
a factor of k g in the hybrid wave function) changes the hybrid states we 
consider to the J PC = 1 h exotic hybrid mesons, those of most interest 
experimentally. In such cases, the decay amplitudes calculated in a LF 
scheme should exhibit the required angular dependence in the CMF of an 
exotic hybrid. In calculations using the standard form of dynamics, one 
should amake sure that the boost symmetry is respected. However, already 
on the basis of our analysis of the non-exotic ++ hybrid decays, we suggest 
that no matter what scheme one uses, a complete set of constraints of special 
relativity should be seriously taken into account in searches for a suitable 
constituent picture. 
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Appendix A 

Basic definitions 



Spinors we use are defined as 



u mp \ = B(p,m) u ox 



v mp x = B(p,m)v ox 



(A.l) 



where the operator B(p,m) 



B(p, m) 



1 



[A + p+ + A_(m + crV 1 )] 



(A.2) 




2S 



represents a boost that changes the mass m at rest into the four-momentum 
p. Spinors of fermions at rest uq\ are 



Mot = v2m 



X+ 




<2m 



and for anti-fermions at rest vq\ are 



'2m 



X- 





-x+ 



(A.3) 



(A.4) 



where x± are two-component spinors, x+ = (1)0) an d X- = (0, 1) . We use 
the convention A± = ^jo^, where ^ = 7 ± 7 3 . 

The integration measure over momenta in a meson, denoted by [k\k2\ or 
just [12], is 



2(2tt) 3 ^ 2(2vr) 3 ^ 2(27r) 3 a;i 2 (l - x 12 ) 2(2vr) 3 P 1 t 



(A.5) 



In terms of the three- vector £12, the integration measure for two constituents 
with the same mass m g , is given by 



dx\od 2 k 



12 a rc 12 



4d 3 £ 



12 



Therefore, 



Similarly, 



2(27r) 3 x 12 (l-xi 2 ) 2(2vr) 3 Mi2 
N pm (kn) = 



IM 



12 



2m n 



N ps (k 12 ) 



^Tr^(l,2)5 p (l,2) 



(A.6) 



(A.7) 



(A.8) 



In the normalization equation for a hybrid meson, we have a three- 
particle integration measure 



n 



dktd 2 k- 



M 3 k n 



d%M qg 



1 167r3k t 2(2tt) 3 M 9 2 (2vr) 3 ^m2 + k 2 ^ ' M\ + kj 



A.9) 



Therefore, 
N hs {k q ,k g ) 



1 Mq 
2m n 



<2m + m 



'/// 



\ 



m 



\ 



2771,0 



(A.10) 
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and 



N hs (k q ,k g ) 



\ 



£Tr s*(l,2,3)Si(l,2,3) 

pol 



(A.ll) 



where X^poi means a sum over two transverse polarizations of a gluon. 





Fig. 11. Hybrid meson decay amplitude into two non-exotic mesons p (the light 
one) and b (the heavy one). 

The decay amplitude of a hybrid into two mesons is 



A(p,b,h) = (-1) 



2 1 g\ f dxucpKj^ f dx^2d 2 K^2 

371(16^)27 X U {1 - X 14 ) J S52(l - ^52) 



x N p N b N h r p (h ±W b (5, 2) [21(1, 2, 3, 4, 5) + 58(1, 2, 3, 4, 5)] 
x [21(1,2,4,5) +55(1,2,4,5)] , 



where 



21(1,2,4,5) = —7^(1, 2, 3, 4, 5)^4(1,2,3,4,5) 

X3 I k 3 =k 4 +k 5 



55(1,2,4,5) = -T B (1, 2, 3, 4, 5)5(1,2, 3, 4, 5) 

%3 



^3=^1+^2 



(A.12) 
(A.13) 



A(l,2,3,4,5) =V^(l,2,3)/ A (Mi 5 ) 
5(1,2,3,4,5) = ^(5, 4, 3)/ A (M? 2 ) 



(A.14) 
(A.15) 



.30 



and the spin factors in the decay amplitude are 

T A (1, 2, 3, 4, 5) =lv[sJ(l,4)5 h (l J 2,3)5 6 t (5,2)5Q C£ ,(5 J 4,3)] , (A.16) 
T B {1, 2, 3, 4, 5) =Ty[4(1,4)S qcd (1, 2, 3)^(5, 2)S h {5, 4, 3)] . (A.17) 

Parts A and B refer to the two arrangements of quarks shown in Fig. ^2 
Sqcd is the spin factor coming from the interaction Hamiltonian of QCD: 



(A.18) 



In both parts of the amplitude, A and B, one has x p = p + /h + , x b 
b + /h + = 1 — x p . In meson p, one has ku = k p , so that M.\i = .M p 



2w + k p 2 , and the following relations hold: 



x u = (\/m 2 + kp 2 + kp)/M p , 



X\ — X\4Xp, 



x u p 



k\ — x\^p + fc„ , 



Z4 = (1 - xu)xp, k£ = (1 - xi 4 )p + , k i =(l-xu)p -k 



(A.19) 
(A.20) 



In meson b, one has £52 = k b , so that A^j, = .M52 = 2ym 2 + k 2 , and the 
analogous relations are: 



0052 = {yrn 2 q + k b 2 + k 3 b )/M b , 
k£=x 52 b + , x 5 = x 52 x b , foj- = x^b 1 - + ki 1 , 

k 2 = (1 - x 52 )b + , x 2 = (1 - x 52 )x b , k 2 = (1 - x 52 )b x - k b . 



(A.21) 
(A.22) 



Evaluating the quarks invariant masses in the hybrid and in the decay 
vertex, one obtains 



Ml 2 = (x 1 + x 2 ) 



kj- 2 + m 2 k 2 + m 2 



xi 



x 2 



(k\ + k 2 ) 2 , 



(A.23) 



and 



M 2 U 



(x 5 + Xi) 



k^ 2 +m 2 

%5 



+ 



k^ 2 + m 2 

Xi 



(k 5 + k 4 



(A.24) 
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The three-vectors: k± 2 in mesons, and kh q and kh g in the hybrid, are 
defined using 



^k 2 + m\ + \Jk 2 + m2 



2 Jl _,_ ^,2 ^2 _|_ ™2 
+ 



«- + mf , ^ + m S _ jw2j (A 25a) 



x 



1 — X 



k = k . 



For mi = m2 = m, one has 



4 (P + m 2 )=4^ ! 4 = -M 2 , 
/ x(l — x) 



k L = K. 



In the part A of the decay amplitude, one has 

kn = M\ 2 /A - ml , 



and 



k„ 



M'( 23 - (M12 + m g f M\ 23 - (Ma - m 



AM 2 



(A.25b) 

(A.26a) 
(A.26b) 

(A.27) 
(A.28) 



123 



where 
M\ 23 = 



ki 2 + m 2 q k^ + m 2 (k 1 + k 2 ) ±2 + m 2 g 
- 1 ' - (fci + k 2 + K3) 



+ 



+ 



Xl 



X2 



\ — X\ — X 2 



Similarly, in the part B, one has 



and 



k„ 



k 2 = M\jA-ml 



M 2 543 - (A4 54 + m g ) 2 M\ A3 - (A4 54 - m g 



AM 2 



(A.29) 



(A.30) 



(A.31) 



543 



where 
Ml 3 = 



£5 



£4 



1 — X5 — X4 



(A.32) 
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Appendix B 

Cross-checks 

All our calculations of spin factors were done using two independent 
methods. One method was to first reduce the spin factors to 2 x 2 matrices 
sandwiched between two component spinors Xa, using Eq. l|T^jl for the gluon 
polarization four-vector e M in A + = gauge. The spin factors were obtained 
from the trace of the product of 2 x 2 matrices. The other method made use of 
the expressions for ^ uu, ^ vv and ^ e*^e v ' . The spin factors were obtained 
using properties of traces of products of 7 matrices. In the evaluation of spin 
factors, we assume that the gluon is massless and has only two degrees of 
freedom. But the gluon acquires an effective mass dynamically. Therefore, 
we used kg = m 2 in the momentum-dependent factor of the wave function. 
In k% and P^ 2 3 m the lattice-inspired spin factors, we checked what happens 
in both cases, i.e., when one inserts fc| = m 2 g or fcf = 0. 

The six-dimensional integrals were carried out using Monte Carlo inte- 
gration (using the procedure VEGAS [231 )• The accuracy of the results of 
integration (standard deviation output from VEGAS) is shown as error bars 
in plots, unless the error is smaller than the size of a point on a plot. The 
VEGAS calculations using C were checked against iterative Gauss quadra- 
ture, also in C, and against a separate FORTRAN program performing the 
same calculations in several representative (but not all regular) cases. 

Appendix C 

Illustrative examples of minimization 

This Appendix provides examples of numerical evidence that we have 
gathered in all cases (many more than given here) we studied. The figures 
show how stddev (standard deviation) or maxdev (maximal deviation), both 
in ratio to the amplitude averaged over the angle 9, change around a global 
minimum when one changes just one parameter in the wave functions or the 
effective Hamiltonian width A. The varied parameter is on the horizontal 
axis, and the deviation on the vertical axis, stddev marked on the right-hand 
scale (plotted with black squares) and the maxdev marked on the left-hand 
scale (plotted with circles). 

The examples demonstrate the dominant feature that the parameters 
fthq and A are strongly correlated with each other and both much larger 
than all other parameters. The first example with the hybrid wave function 
with the spin factor u^v and N ^ 1. The example illustrates our argument 
for that the parameters fih q and A must be both much larger than all other 
parameters: the minimum resembles one side of a broad valley open toward 
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large values. The remaining four examples concern the cases with various 
spin factors and always N = 1. The examples illustrate that j3h q and A are 
both strongly correlated and much larger than all other parameters. 
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Fig. 12. Variation of rotational symmetry violation versus changes of parameters 
in the wave functions in the case a) in Fig. El i.e. Sh — u^v, with extra factors 
N 5^ 1, decay into two J PC = ++ mesons. The optimal values of the parameters 
are given in the first column in Table El The arrow marked "ASA" points toward 
the optimal value of a parameter. The last plot shows the amplitude itself for the 
optimal parameters. 
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Fig. 13. Variation of rotational symmetry violation versus changes of parameters 
in the wave functions in the case b) in Fig. El i.e. Sh = u^v, with N — 1, decay into 
two J PC = ++ mesons. The optimal values of the parameters are given in the 
second column in Table 03 The arrow marked "ASA" points toward the optimal 
value of a parameter. The last plot shows the amplitude itself for the optimal 
parameters. 
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Fig. 14. Variation of rotational symmetry violation versus changes of parameters 
in the wave functions in the case b) in Fig. i.e. Sh — u$v, with N = 1, decay 
into two J PC = 0~ + mesons. The optimal values of the parameters are given in 
the second column in Tabled The arrow marked "ASA" shows value of parameter 
found in minimization. The last plot shows the amplitude itself for the optimal 
parameters. 
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Fig. 15. Variation of rotational symmetry violation versus changes of parameters 
in the wave functions in the case f) in Fig.0 i.e., Sh = uGPv, with N = 1, decay 
into two J PC = ++ mesons. The optimal values of the parameters are given in 
the sixth column in Table 0] The arrow marked "ASA" points toward the optimal 
value of a parameter. The last plot shows the amplitude itself for the optimal 
parameters. 
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Fig. 16. Variation of rotational symmetry violation versus changes of parameters in 
the wave functions in the case f) in Fig.EHS'/j = uGPv, with N = 1, decay into two 
" + mesons. The optimal values of the parameters are given in the sixth 
column in Table El The arrow marked "ASA" points toward the optimal value of a 
parameter. The last plot shows the amplitude itself for the optimal parameters. 



